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Abstract 

A non-equilibrium theory of isothermal and diffusionless evolution of incoherent interfaces 
within a plastically deforming solid is developed. The irreversible dynamics of the interface are 
driven by its normal motion, incoherency (slip and misorientation), and an intrinsic plastic flow; 
and purely by plastic deformation in the bulk away from the interface. Using the continuum 
theory for defect distribution (in bulk and over the interface) we formulate a general kinematical 
framework, derive relevant balance laws and jump conditions, and prescribe a thermodynami- 
cally consistent constitutive/kinetic structure for interface evolution. 

keywords: Continuous distribution of dislocations, Finite strain elasto-plasticity, Incoherent interfaces, 
Interface evolution. 

1 Introduction 

The motivation for the current work is derived from the processes in material evolution where a 
moving interface plastically deforms the bulk material region, for example recrystallization and 
impact induced plasticity [101 125] 134] . The interface is taken to be a sharp surface separating 
two distinct regions such as different phases (during phase transition), different crystals (in poly- 
crystalline materials), and differently oriented single crystals (grain boundary), or a wave front 
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during dynamic deformation. Even after we assume the processes to be isothermal and diffusion- 
less and make simplifying assumptions about the bulk and the interface, the rich dynamics of an 
interface offers a wide gamut of problems to the material scientist |43] and the mathematician [19] 
alike. The challenge on one hand is to construct physical models which are amenable to experi- 
mental verification and numerical implementation, and on the other hand to analyze the resulting 
partial differential equations for their well-posedness and properties of the solutions. 

The structural nature of the interface is characterized on the basis of its behavior upon relaxation 
of local stresses. We call an interface incoherent if, upon relaxation, it is locally mapped into two 
disjoint configurations. Otherwise, we call it coherent. An incoherent interface in an otherwise 
defect-free solid, after stress relaxation, will result into two separate solids [H[33]. Incoherency is 
expressed in terms of the incompatibility of the distortion field and leads to interface dislocation 
density as a smeared-out defect distribution (cf. Bilby and coworkers [21 El E] and Ch. 2 of [43]). 
The interfacial dislocation density along with its bulk counterpart contributes to the Burgers vector 
for arbitrary closed curves (crossing the interface) in the body. If the interface is coherent then its 
defect distribution, and consequently its contribution to the net Burgers vector, vanishes identically. 

We consider plasticity to be a purely dissipative phenomenon driven by irreversible changes in 
the microstructure. Even in the bulk, away from the interface, the evolution of plastic flow is a 
complicated non-linear problem coupled with elasticity and non-local microstructural interactions. 
Many of the underlying mechanisms remain poorly understood [35] and it is becoming increasingly 
necessary to develop the theory at a microstructural level. One related concern is to understand the 
plastic behavior at internal boundaries separating different phases or grains [27]. The plastic flow 
behavior at such boundaries will depend on interface motion, relative distortion of the neighboring 
grains, and the local shape (for example orientation and curvature). It is clear that a theory for 
plastic flow at the interface cannot be, in general, modeled along the same lines as the theory 
associated with the bulk. 

Interfaces in solids, with an associated energy density, have been well studied in the context of 
continuum thermodynamics. We note, in particular, the earlier work done to obtain equilibrium 
conditions for coherent /incoherent interfaces within elastically deforming solids [26} [3TJ l32j 133] . 
These conditions were obtained by minimizing the total energy (bulk and interfacial) under ap- 
propriate variations in the domain (see Remark 12.31 for further discussion). Gurtin and coworkers 
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[6j [181 [SOI [22] extended these results to dynamic interfaces and demonstrated the validity of equilib- 
rium interfacial conditions in wider settings than were previously considered. Their methodology 
relies upon a version of the virtual work principle where contributions from configurational forces 
were considered in addition to those from classical forces. All these theories, however, assume 
the bulk surrounding the interface to be defect-free and thus neglect any possibility of interaction 
between interfacial and bulk defect densities. They therefore fall short of modeling the behavior 
of interfaces in a plastically deforming medium. On the other hand, some recent strain gradient 
plasticity models with interface energies dependent on (infinitesimal) plastic strains [1] [H] [T5] incor- 
porate interfacial flow rules along the same principles as those in the bulk. These relations furnish 
boundary data for the bulk equations. While restricting themselves to infinitesimal strains, these 
models also neglect any coupling with other processes (for example the motion and the relative 
distortion of the interface) . 

Our aim is to generalize the above mentioned works by developing a continuum theory for 
interface evolution in a plastically deforming solid under isothermal and diffusionless conditions. 
Both the bulk and the interface are assumed to possess a continuous distribution of defects, whose 
density is related to the local elastic and plastic distortion maps. The role of a relaxed manifold 
is emphasized in the multiplicative decomposition of deformation gradient in the bulk and at the 
interface. We restrict our developments to the point of positing specific kinetic laws and therefore 
stop short of formulating complete boundary-initial-value problems. We however provide a detailed 
description of the associated kinematics, derive all the necessary balance laws and jump conditions, 
and use physical and material symmetries to restrict the form of constitutive/kinetic relations. 
In particular, we derive local dissipation inequalities and highlight the interplay between various 
dissipative mechanisms and the associated driving forces. The bulk behavior in this paper is 
modeled after our recent work |16t ITT] on bulk plasticity. 

The central results in this paper are: 

(i) The multiplicative decomposition of the interface deformation gradient is equivalently given 
in terms of two sets of (interfacial) elastic and plastic distortions, cf. (|2.87p . Both of these coincide 
for coherent interfaces. 

(ii) The relation between an incoherency tensor and true interface dislocation densities, cf. 
(1332]) and (f3T33D . 
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(Hi) The relationship between bulk and interface dislocation densities given in (|3.44|) . which 
also highlights the fact that interface dislocation density, unlike the bulk dislocation density, does 
not have a vanishing divergence. 

(iv) The dissipation inequality (|4.36p arising due to interface motion, plastic flow at the interface, 
and change in relative distortion across the interface. This inequality demonstrates the underlying 
coupling between the interface motion, the tangential plastic distortion of the neighboring grains, 
and the relative tangential distortion of the grains. It provides a starting point for developing 
kinetic laws governing the out of equilibrium thermodynamic process. Otherwise, in thermodynamic 
equilibrium, it furnishes additional balance laws to be satisfied at the interface. 

(v) The restrictions on the form of kinetic laws, (|4.62p - (l4.64h . due to various symmetries in the 
model. 

Our work furnishes the pre-requisite information about kinematics, dissipation, and the basic 
requirements for constitutive equations needed for the formulation of complete boundary-initial- 
value problems in the study of dynamic incoherent interfaces within plastically deforming solids. 

We have divided this work into three parts. In the first, we prepare the necessary background 
for studying the thermodynamics of energetic interfaces within a bulk medium. The second part 
is concerned with the interface dislocation density as a measure of defect distribution over the 
interface and its relation with the bulk dislocation density. The final part deals with the energetics 
and kinetics of incoherent interfaces moving within plastically deforming solids. We make certain 
constitutive assumptions about the nature of interfacial energies and use them to evaluate the 
net dissipation at the interface. Motivated by the dissipation inequality, and exploiting various 
symmetries of the physical space and the material, we formulate restrictions on kinetic laws at the 
interface. 

2 Preliminaries for the theory 

In the following we prepare the ground work for the next two sections. Our discussion on the 
kinematics and thermodynamics of surfaces, in Subsections 12.21 and 12.31 is largely based upon 
the work of Gurtin and coworkers [U [T8| 1201 [22] and Silhavy |39[ |4"0| |4"T] . Our derivation of the 
interface dissipation inequality, cf. (|2.70p or (I2.76p . is however different and appears to be new. 
Similar relations were obtained in [61 [T51 [22] within the framework of configurational mechanics 
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2.1 Three-dimensional continuum 

The translation space of a real three-dimensional Euclidean point space £ is denoted by V. Let Lin 
be the space of linear transformations from V to V (second order tensors) . The groups of invertible 
tensors, orthogonal tensors, and rotations are denoted by InvLin, Orth, and Orth + , respectively. 
The spaces of symmetric, symmetric positive definite, and skew tensors are represented by Sym, 
Sym + , and Skw, respectively. The determinant and the cofactor of A G Lin are denoted by J a arid 
A*, respectively, where A* = JaA.~ t if A € InvLin (superscripts T and —1 denote the transpose 
and the inverse, respectively, and A~ T = (A~ l ) T ). The space Lin is equipped with the Euclidean 
inner product and norm defined by A • B = tr(AB T ) (B G Lin) and | A| 2 = A • A, respectively, 
where tr(-) is the trace operator. 

We use both indicial as well as bold notation to represent vector and tensor fields. The compo- 
nents in the indicial notation are written with respect to the three-dimensional Cartesian coordinate 
system. Indices are denoted with roman alphabets appearing as subscripts. Summation is assumed 
for repeated indices unless stated otherwise. Let be the three-dimensional permutation symbol; 
it is 1 if {i,j, k] is an even permutation of {1,2,3}, —1 if it is an odd permutation, and if any 
index is repeated. 

Let n r C £ and Kt C £ be the reference configuration and the spatial (or current) configuration 
with translation spaces V Kr C V and V Kt C V, respectively. There exists a bijective map x between 
K r and Kt; therefore for every X € K r and time t we have a unique x € Kt given by 

x = *(X,t). (2.1) 

We assume x to be continuous but piecewise differentiable over K r and continuously differentiable 
with respect to t. 

The derivative of a scalar valued differentiable function of tensors G : Lin — > R (where R is the 
set of all real numbers) is a tensor Ga defined by 

G(A + B) = G(A) + Ga ■ B + o(|B|), (2.2) 

where [gj — > as |B| — > 0. Similar definitions can be made for vector and tensor valued 
differentiable functions (of scalars, vectors, and tensors). In particular, if the domain of a function 
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is K r we denote the derivative by V; and if it is Kt then we write grad for the derivative. Such 
functions are called fields. The divergence and the curl of fields, on n r , are defined by (for w € V) 

Divw = tr(Vw), (Curlw) • c = Div(w x c), (2.3) 
(Div A) • c = Div(A T c), and (Curl A)c = Curl(A T c) (2.4) 

for any fixed c G V. Similar definitions hold for fields on K t ; in this case we denote divergence and 
curl by div and curl, respectively. The material time derivative is the derivative of a function with 
respect to time for fixed X; we denote it by a superimposed dot. 

The particle velocity v € V Kt is defined as v(X, t) = x(X.,t). If x is differentiable at X, then 
the deformation gradient F € InvLin : V Kr —> V Kt exists at X and is given by F = V%. We assume 
v and F to be piecewise continuously differentiable over n r ; they (and their derivatives) are allowed 
to be discontinuous only across the singular surface. 

2.2 Singular surface 

A singular surface (or interface) is a two dimensional manifold in the interior of K r (or Kt) across 
which various fields (and their derivatives) may be discontinuous, which otherwise are continuous 
in the body. A singular surface in n r is given by 

S r = {X € K r : 0(X, t) = 0}, (2.5) 

where ^ is a continuously differentiable function. The unit normal to the surface and the normal 
velocity are defined by 

t/(X, t ) = - ] A f , (2.6) 

respectively; the derivatives being evaluated at the surface. The projection tensor 1(X, t) which 
map vectors in V Kr to vectors in T$ r (X) C V Kr , where T$ r (X) is the tangent space at X G S r such 
that N _L T Sr (X), is given by 

1 = 1 -N<g>N, (2.7) 
where 1 is the identity tensor in Lin. Note that t T = 1 and 11 = 1. 
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The jump in a discontinuous field (say VP) is denned on the singular surface and is denoted by 

m = y+ - (2.8) 

where ^ + and are the limit values of $f as one approaches the singular surface from either side. 
The + side is the one into which the normal to the surface points. Let <3? be another piecewise 
continuous field. The following relation can be verified by direct substitution using (j2.8|) : 

[*¥J = [«&!(*> + ($)[fl, (2.9) 

where 

<*> = 1 • (2-10) 

Derivatives on the surface We first introduce the general idea of derivatives on manifolds 
embedded in a higher dimensional space (see for example \39\ I41j). Let Ai C Lin be a manifold in 
the space of tensors. The derivative of a scalar valued differentiable function of tensors g : A4 — > R 
is a tensor g A defined by (for {A,B} € M) 

g(A + B) = g(A) + g A ■ B + o(|B|), such that g A P(A) = g A , (2.11) 

where -^j]fp — > as |B| — > 0, and P(A) is the orthogonal projection onto the tangent space of M. 
at A. Similar definitions can be made for vector and tensor valued functions on manifolds. 

Let f, v, and A denote a scalar, vector, and tensor valued field, respectively, on S r . They are 
differentiable at X € S r if they have extensions /, v, and A to a neighborhood of X in K r which 
are differentiable at X. The surface gradients of f , v, and A at X G S r are defined by 

V s f(X) = 1(X)V/(X), (2.12) 
V 5 w(X) = Vv(X)l(X), and (2.13) 
V 5 A(X) = VA(X)1(X). (2.14) 

In the rest of the paper we will use the same symbol for both the surface field and its extension. 
We define the surface divergence of v as a scalar field Div 5 w; and of A as a vector field Div s A 
given by 

Div 5 v = tr(V 5 v) and 

c • Div 5 A = Div s (A T c) (2.15) 



for a fixed c £ V. Moreover, we call v (or A) tangential if lv = v (1A = A) and A superficial if 
Al = A. 

Define the curvature tensor L by 

L = -V 5 N. (2.16) 

It is straightforward to verify that L = ~L T (use (|2.6P i ) and LN = 0. Therefore, N is an eigenvector 
of L with zero eigenvalue. Since L is symmetric, the spectral theorem implies that it has three real 
eigenvalues with mutually orthogonal eigenvectors. Let the two nontrivial eigenvalues be d and 
£2 with eigenvectors in T$ r ■ The mean and the Gaussian curvature associated with the surface are 
defined as 

# = ^(Ci + C 2 )andK = CiC2, (2.17) 

respectively. 

A function Lp : (t — e, t + e) — > n r , e > is said to be a normal curve through X £ S r at time t 
if for each r £ (t — e, t + e), (p(r) £ S r and 

<p'(t) = U(cp(t),tM<p(t),t), (2.18) 

where the superscript prime denotes the derivative with respect to the scalar argument. Define the 
normal time derivative of a field on S r by (cf. §11.3 of |44| and §179 of |45] ) 

HX , t) = M*pA . (2 . 19) 

dT r=t 

It represents the rate of change in v with respect to an observer sitting on S r and moving with 
the normal velocity C/N of the interface. As an example, on differentiating (|2.6P i and using the 
definitions for surface divergence and normal time derivative, we obtain 

N = -V S U. (2.20) 

Therefore, evolving surfaces S r are parallel if and only if U is constant over S r at any fixed time. 

Compatibility conditions The continuity of deformation field x(X, t) across S r furnishes the 
following jump conditions for the deformation gradient and the velocity field (cf. Ch. II of [H] and 
Ch. C of 05]): 
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[F] = k ® N and (2.21) 
[v]+[/[F]N = OVXeS r , (2.22) 

where k £ V Kt is arbitrary. For U ^ these relations can be combined to eliminate k, 

l7[Fj = -[v]®N. (2.23) 

Singular surface in the current configuration The image of the singular surface S r in the 
current configuration is given by 

s t = {x € Kt : </>(x, t) = 0, with ^(*(X, t),t) = 0(X, t)}. (2.24) 

The scalar function ^ is continuous but, in general, only piecewise differentiable with respect to its 
arguments. The derivatives of tp can suffer jump discontinuities at st- Differentiate ip(x(X-,t),t) = 
0(X, t) (away from st) with respect to X (at fixed t) and t (at fixed X), and then restrict the result 
to the surface, to obtain respectively, 

= (F ± ) T (gradV') ± and 
0=(grad^) ± -v ± + , (2.25) 

where ± indicates that either of + or — limit of the field can be used to satisfy the equation (due 
to smoothness of <fi across the singular surface), and ^ indicates the partial derivative of ip with 
respect to t at fixed x. Substitute (|2.6p into f|2.25|) to get 

(grad^) ± (F ± )" T N (F ± )*N 



and 



(gradV^I |(F±)- T N| |(F±)*N| 

1 /d^N* _ (grad^) ± U 



|(grad^)±| \dtj |(gradV) ± | V + |(F±)- T N|' (2 " 26) 
The compatibility relations (|2.21|) and (|2.22|) yield the + and — value of the expressions on the 
right hand sides above identical, cf. f|2.33j) below. This leads us to define the normal to the surface 
st and the spatial normal velocity by, cf. (|2.6p . 



(2.27) 



= (gradVQ^ 
|(grad r 4>) ± \ 



1 /<9^ xd 



|(grad ip) ± \ V dt 
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respectively; we obtain 

(F ± )" T N (F±)*N 
m ~ |(F±)~ r N| ~ |(F±)*N| 

M = m - V±+ 1(F^N| - (2 ' 28) 
The projection tensor 1(X, t) which map vectors in V Kt to vectors in T St (x) C V Kt , where T St (x) 
is the tangent space at x = %(X, t) such that e _L T St (x), is given by 

l = l-e<g>m. (2.29) 

Surface deformation gradient and normal velocity For a continuous motion across the 
surface, i.e. [%(X, t)J = for X 6 S r , we define the surface deformation gradient F and the surface 
normal velocity v on S r as j22[ [39] 

F = V s x and v = X- (2-30) 

It is then easy to check that 

F = F ± l and v = v ± + C/F ± N. (2.31) 

Tensor F satisfies detF = 0, which can be verified using (I2.3ip i and det 1 = 0. Moreover, we 
have from ([231^ and $2?28b i. 

FN = and F T ii = 0. (2.32) 

Therefore, Fl = F and IF = F. The cofactor F* of F is defined by F*(a x b) = Fa x Fb for 
arbitrary vectors {a, b} € V Kr . Let {ti,t2} £ 7s r (X) be two unit vectors such that {ti,t2,N} 
forms a positively oriented orthogonal basis at X. Then 

F*N = F*(ti x t 2 ) = Ft! x Ft 2 

= F ± ti x F ± t 2 

= (F ± )*N, (2.33) 

where in the third equality we have used (12.3^ . On the other hand, employ (|2.32j) i to conclude 
that F*t a = (a = 1,2) and hence F*l = 0. Therefore, F* remains non-zero because (F±)*N 
does not vanish. According to (|2.33j) . |F*N| is equal to the ratio j of the infinitesimal areas (on the 
singular surface) in the current and the reference configuration. Use F* = F*l = F*N <8> N to write 

F*=j(n®N). (2.34) 
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Hence |F*| = |F*N| = j. 

Following Penrose [37] we define a unique tensor F _1 , the pseudoinverse (or the generalized 
inverse) of F , such that 

F _1 F = 1 and FF _1 = I, (2.35) 

which also satisfies 

F- 1 = (F ± )" 1 1, (2.36) 

as can be checked by direct substitution. 

For F^ G InvLin there exist unique tensors R 1 * 1 £ Orth + and U 1 * 1 € Sym + such that F^ = 
R ± U ± . For a non-invertible tensor F there exists a unique positive semidefinite tensor U € Sym 
and a (non-unique) orthogonal tensor R G Orth such that F = RU. These statements follow from 
the polar decomposition theorem for invertible and non-invertible tensors. Recall (j2.31[) i to write 
F = R ± U ± 1. Tensor U thus satisfies U 2 = F T F = 1U ±2 1. Define R = FIT 1 , where IT 1 is the 
pseudoinverse of U such that U _1 U = UU^ 1 = 1. Tensor M. is unique and satisfies 

M. T R = 1 and MM T = 1. (2.37) 

Moreover, tensor R in the polar decomposition for F is related to M as Rl = R. The expression 

F = MJU (2.38) 

provides a decomposition for F into unique tensors. 

The surface gradient of normal velocity can be calculated from (|2.31() -? 

V s v = F ± l -F ± N(g)N- C/F ± L, (2.39) 

where, in addition to the definitions of surface gradient and normal time derivative, we have used 
(235} and Vv ± = F ± . Employ fl23Di and 

i = -N<8>N-N<8>N (2.40) 

to rewrite (|2.39p as 

V s v = Fl - UFL. (2.41) 
Consequently it is only for a flat interface (L = 0) that we have V^v = Fl (compare with Vv = F). 
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Remark 2.1. Let g be a scalar function on the interface given by g = g(¥, N). The arguments of 
g satisfy FN = and \N\ = 1 and therefore form a submanifold, say G, of Lin x V. The partial 
derivatives <?f and <?n (with respect to F and N, respectively) are evaluated using an extension g of 
g and restricting the result to G. Extension of g is any smooth function defined over Lin x V such 
that it is equal to 5 on G. These partial derivatives lie in the tangent space of G and hence satisfy 
(cf. (12. 1 1 f) : for a proof see Appendix B of [41 J ) 

g ¥ N + F<? N = and N • g N = 0. (2.42) 

In the rest of the paper we will use same notation for the function and its extension. 
Remark 2.2. (Derivative of j) Use (|2.33p and (|2.34|) to obtain 

j 2 = (Fti x Ft 2 ) • (Fti x Ft 2 ), (2.43) 

where ti and ti 2 are functions of only N; i.e., they are arbitrary orthonormal vectors orthogonal to 
N. To find partial derivative j'n fix F in (|2.43p and differentiate it on a one-parameter curve in the 
space of all unit vectors satisfying FN = 0. Apply the definition of cofactor and use (I2.32p i to get 
j N = 0. Therefore, by (ET421) . j ¥ N = 0. On the other hand, differentiating j = |(F ± )*N| for fix N 
yields 

jw = JF" T . (2.44) 
Hence the normal time derivative of j is given by (compare with Jp = J^FF -1 • 1) 

3'=jFF _1 -I. (2.45) 
2.3 Balance laws and dissipation inequality 

Assuming a purely mechanical environment and isothermal heat flow we obtain balance laws for 
mass and momentum, and the dissipation inequalities both for material points on the interface and 
away from it. We do not state the balance of energy since it is used, under isothermal conditions, 
only to calculate the net heat flux during the dissipative process. 

Surface divergence theorem and surface transport theorem In addition to divergence 
and transport theorems for piecewise smooth fields on K r (see for example Ch. 3 of [38]) we 
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will repeatedly use the following theorems for fields defined on S r . For a vector field w € V Kr 
continuously differentiable on S C S r 

/ w • vdL = / (Div 5 w + 2ffw • N)dA, (2.46) 
JdS Js 

where v € T$ r is the outer unit normal to the closed curve dS bounding S such that {N, t} forms 
a positively oriented orthonormal basis on dS with t as the tangent vector along dS. Moreover if 
w is tangential, i.e. Iw = w, then w • N = and f)2.46j) reduces to 

w • udL = i Div 5 wdA. (2.47) 
as Js 

The surface transport theorem for an evolving surface S within a fixed region f2 such that 
dS C dQ is given by [23] 



d_ 

If 



/ wdA= / (w - 2UHw)dA - / wll cot 6dL, (2.48) 
Js Js JdS 

where 6 = arccos(N • N) and N is the outward unit normal on dQ. If 0, is arbitrary then we can 

always choose with d£l such that N • N = at all X £ dS i.e., orient d£l in such a way that it is 

orthogonal to S at all points on dS (cf . Figure [2]) . With this choice (|2.48H reduces to 



d 
If 



[ wdA = f (w - 2UHw)dA. (2.49) 
Js Js 



Similar theorems hold for scalar and tensor fields on S r . 

Conservation of mass Assume no net mass transfer in an arbitrary volume of K r . Also assume 
that there is no additional mass density associated with S r . The statement of conservation of mass 
then reduces to [35] 

p r = VX G K r \S r , (2.50) 
where p r is the referential mass density of the bulk, and 

U{p r } = VX € S r (2.51) 

i.e., either the normal velocity vanishes or the referential mass density is continuous 
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PN 



Figure 1: Various forces acting on an arbitrary volume f2 C n r . 



Balance of momentum The balance laws for linear and angular momentum can be either stated 
as Euler's postulates or can be deduced from the first law of thermodynamics [38]. Let Q be a three- 
dimensional open subset of n r with boundary dtl such that S = f2 n S r is nonempty and dS C d£l. 
Let N E V Kr and v E T$ r be unit vectors normal to d£l and dS, respectively. Let P E Lin be the 
bulk Piola stress and b E V the specific body force vector. We assume the existence of a contact 
force between two subsets of S r along the curve of contact, which can be expressed in terms of 
a linear map (given by interface Piola stress P E Lin) acting on the normal to the contact curve 
|21j . If there are no body forces associated with the singular surface then the balance of linear 
momentum for Q is given by (see Figure [1] where all the forces are shown) 



Let P be superficial, i.e. PI = P. This is motivated from the last term of the above equation 
where PN <g> N does not contribute to the net force (since N is orthogonal to u) and therefore can 
be assumed to vanish without loss of generality. Integral equation (|2.52j) can be localized, using 
the transport and divergence theorems, to |18[ [22] 




(2.52) 



p r v = Div P + p r h VX £Kr\S r and 
Uprfv] + [PJN + Div 5 P = VX E 5, 



(2.53) 



(2.54) 



where we have also used (|2.50p and (|2,51|) , 
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The balance of angular momentum is given by 

/ p r r x vdV = r x PNcL4 + / Pr v x bdV + / r x PudL, (2.55) 
Jn Jan Jn Jas 



d_ 

di 



where r = x — xo and xo € £ is arbitrary. On using transport and divergence theorems and the 
equations of balance of mass and linear momentum, it localizes to [18\ [22] 

PF T = Fp T vx e Kr y Sr and ^ 5g ^ 

PF T = FP T VX £ S r . (2.57) 

Equation (f23D) implies that IP = P. Indeed, use IF = F to get 1PF T = PF T . The desired result 
follows upon using ()2.31j) i , the invertibility of F 1 * 1 , and PI = P. 

The interface Cauchy stress T € Lin is a superficial tensor (Tl = T) which satisfies 



Yvdl = / FvdL (2.58) 

It JL t 

for Lt (with normal u) and It = x(Lt) (with normal u) as curves on the referential and spatial 
singular surface, respectively. Let {u, t, N} be a positively oriented orthonormal basis on S r . Define 
t £ T St by tdl = F^t. The triad {P,t,m} then forms a positively oriented orthonormal basis on 
St, where m is given by (12 .28[) i . Hence vdl = j~ 1 (F :il t x (F ± )*N)(iL, which on repeated use of the 
definition of cofactor simplifies to 

vdl = j¥~ T udL. (2.59) 

Stresses P and T are therefore related as (compare with P = J^TF _T , where T S Lin is the bulk 
Cauchy Stress) 

P = jT¥~ T . (2.60) 

The balance laws in the spatial configuration, equivalent to (|2.53p . (|2.54p . (|2.56p . and (]2.57p . 
are given by 

pv = div T + pb, T = T T Vx € K t \ s t , (2.61) 
j s lvj + pTjn + div 5 T = 0, and T = T T Vx G s t , (2.62) 

where p is the mass density with respect to Kt and j s = p r Uj~ l . 
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Dissipation inequalities Let ^ and <£> be the free energy densities per unit volume of K r and 
per unit area of S r , respectively. Assume that S r has zero body force and kinetic energy density 
For an arbitrary volume f2 C K r , with S = n S r nonempty and dS C d£l, the mechanical version 
of second law of thermodynamics (under isothermal conditions) yields 

/ Pr b-vdV+ [ PN -vdA + [ Fv-vdL-— I ( * + \p r |v| 2 ) dV - ^ / <S>dA > 0. (2.63) 
Jn J do, J as dt J n \ 2 J dt J s 

A comment is in order for the term representing the power due to interfacial stress. At every 

point on the curve dS the contact force (between the surfaces divided by the curve) is given 

by Pi/ and the rate of change in displacement, with respect to an observer sitting on S (at the 

considered point) and moving with velocity UN, is given by v. The change in displacement apparent 

to the observer sitting on S but moving tangentially to the interface will depend on the chosen 

parametrization and so will the resulting power. This is undesirable and therefore we use only w 

to calculate the power expended at the interface. Gurtin and coauthors [U 1201 [22] have imposed 

invariance with respect to tangential velocities in their formulation of configurational balance laws. 

This is equivalent to the requirement of invariance under re-parameterizations of the interface. 

Our viewpoint is different: We require (a priori) the mechanical power balance to be invariant 

under re-parametrization and write it in a form that satisfies this invariance automatically. Thus 

this requirement is automatically satisfied in the present formulation and accordingly yields no 

non-trivial information. 

Before we proceed let us clarify the nature of interfacial stresses. The interface stress P, in 
contrast to the bulk stress, does not act on a fixed set of material points but rather on material 
points momentarily occupying the surface S. This is in accord with the mechanism responsible for 
surface tension in liquids. As the surface area increases, interstices are generated which are filled by 
molecules from the bulk liquid. In this way the surface tension remains sensibly constant while the 
surface area expands. Thus the matter occupying the surface does not actually stretch. Instead, 
the surface changes its area due to the continuous addition of mass. This physical situation stands 
in contrast to the treatment of surface tension in conventional continuum mechanics, in which the 
surface is regarded as a material surface if the motion of the liquid, regarded as a closed set, is 
continuous. In the conventional interpretation, surface tension is then a conventional force system 
acting on a persistent set of material points. However, in the actual physical situation, the surface is 
not material in the usual sense. Our framework accommodates such mechanisms while retaining the 
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conventional interpretation of force. The contribution to mechanical power from interface stresses 
(as in ()2.63|) above) is consequently obtained not by its action on material velocities but on v. 
Using the transport, divergence, and localization theorems, (12.63H reduces to 



P • F - \P > VX € K r \ S r and (2.64) 

+ Div 5 P T w - ($ - 2UH$>) > VX € S r , (2.65) 



PN • v + U \V + ip r |v| 2 

where, in obtaining (|2.64|) . we have also used balances of mass and momentum. We now rewrite 
(|2.65p using the identities 

[PN ■ vj = [PJN • v - U{PN ■ FN] , (2.66) 

[v • vj = 2[v] • v + C/ 2 [|FN| 2 ], and (2.67) 

Div s P T v = Div 5 P • v + P • V 5 v. (2.68) 

Here ([2T66]) and (|2U7]) can be verified with the help of §2J)§, (pT22]) . and (pl3"T]) o while (|2"U8l) follows 
from the chain rule of differentiation. These identities, in addition to (|2.5ip and (|2.54p . reduce 
([235]) to 

U (n ■ [E]N + ^V[|FN| 2 jJ + P • V s v - ($ - 2C/#$) >0VXeS r (2.69) 
or equivalently (on substituting V s w from (|2.4ip and 2H = tr L) 

U (n • [EJN + ^ 2 p4|FN| 2 ]^ + C/E • L + P • Fl - $ > VX 6 5 r , (2.70) 

where 

E = *1 - F T P and (2.71) 

E = $1-F T P (2.72) 

are bulk and interface Eshelby tensors defined over K r \S r and S r , respectively. Dissipation inequal- 
ities (|2.64p and (|2.70p . in addition to balance laws for mass and momentum, should be satisfied for 
every process. 

Remark 2.3. The present setting differs from that of Gurtin [6] 122] [20] as we do not consider any 
explicit contribution from configurational forces in the global dissipation inequality (I2.63P (compare 
with Equations (21-6) and (21-19) in [20]): the final results however coincide. We demonstrate this 
by assuming, for now, the interface energy density to be of the form $ = $(F, N). Such energies 
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Figure 2: A portion of d£l intersecting with a portion of S r along dS. N is the outward unit normal 
to d£l and N to 5 r . t is the unit tangent along dS. Since fi is arbitrary we can choose d£l such 
that t is parallel to c and N coincides with v, where c £ T$ r is defined in Remark 12.31 and v £ Tg r 
is the outward normal to dS. This ensures that N • N = and c • v = (cf. paragraphs following 
(gaSD and (j2J5]) ). 



have been well studied in the contexts of phase equilibrium with interfacial energy [22 1 133 1 141]. The 
surface stress is given by P = 5>f1. We can then obtain $ = P - Fl — c • N, where <c = F t $fN — 4>pj 
is tangential. Substituting N from ()2.20p and using the chain rule of differentiation yields 



Fl - U Div 5 <c + Div 5 (C/c). 



(2.73) 



Substituting it in (|2.7Up we get 
1 



[7 ( N • [EJN+ -C/ p r [|FN| J ) +?7E-L + C/DiV b c-Div 5 (C/c) > VX € S r . 



(2.74) 



The term Div (U<c) drops out of the inequality. Indeed after integrating (|2.74j) over S and applying 
surface divergence theorem (|2.47p this term takes the form 



U c • i/dL, 



dS 



(2.75) 



where v € T$ r is the exterior unit normal to dS. Since O is arbitrary and c is tangential, we can 
choose 30 such that c • v = (i.e., orient d£l such that c is parallel to the tangent at every point 
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on dS, cf. Figure [2]). Upon localization of the resulting integral inequality we are finally led to 

U (n • [E]N + ^C/ 2 pr[|FN| 2 ] + E • L + Div 5 cj > VX € S r (2.76) 

as a necessary condition for (|2,63p . The coefficient of U is the net driving force for the motion of a 
coherent interface between two bulks phases. This coincides with the result obtained by Gurtin, cf. 
Equations (21-10a) and (21-26) in [20]. The configurational shear r appearing in those equations 
from [20] is equal to <c (see §2 of [5T] in this regard). At thermodynamic equilibrium the driving 
force vanishes thereby furnishing a balance relation to be satisfied at the interface. Such relations 
were also obtained via energy minimization |26[ \31\ [321 [33] . Ours is a dynamical theory, whereas 
results coming from energy minimization are really only relevant at equilibrium, and even then 
only for stable equilibria. 

2.4 Elastic plastic deformation 

The idea of stress-free local configurations is central to our theory. We assume both the bulk and 
the interface stress to be purely elastic in origin, wherein the deformation is measured with respect 
to the stress-free configuration. In a recent paper |16] we demonstrated, using the mean-stress 
theorem, that it is always possible to obtain a locally stress-free state (under equilibrium and in the 
absence of external forces) by cutting m into parts with arbitrarily small volume. Moreover, if these 
sub-bodies cannot be made congruent in absence of any distortion then they do not form a connected 
set in a Euclidean space. The material is then said to be dislocated with no global differentiable 
map from n t to the disjoint set of sub-bodies [H [281 US EH [36] . The union of these unstressed 
sub-bodies is a three-dimensional non-Euclidean smooth manifold, say A4. A local configuration 
in A4 is identified with the local tangent space, denoted by Ki- The local map from to V Kt is 
represented by H £ Lin. The absence of a global differentiable map renders H incompatible and 
therefore, unlike F, it cannot be written as gradient of a differentiable map. The incompatibility 
of H implies the existence of a continuous distribution of dislocations over Kt (see the next section 
for details). 

The argument used for the existence of stress-free local configurations in [16] assumes smooth- 
ness of bulk stress. If the stress field is non-smooth only over a set of measure zero, the stresses can 
still be relaxed on neighborhoods arbitrarily close to the singular region and therefore everywhere 
except over the set of zero measure. If singular regions have stresses associated with them, for 
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example the surface stresses discussed above, then they also need to be relaxed. In the following 
we show that this can be done under equilibrium and vanishing external forces if the surface st is 
cut into infinitesimal areas. 

To this end consider an arbitrary subsurface s C St and assume T to be continuously differen- 
tiable over sj. A simple calculation (using (|2.62p i without the inertial term) then yields 

f = - / Tda = — I [ p(x<g> [Tjn+ [T]in<g> x)da + / p(x ® Ti> + Tu <g> x)dl 1 , (2.77) 
a Js 2o [Js Jds J 

where T is the mean interface Cauchy stress and a is the area of s. The mean stress therefore 
vanishes if there are no external forces on s. According to the mean value theorem, there exists 
x £ s such that T(x, t) = T (= 0). Let the area a become arbitrarily close to zero. Then, by 
continuity of T, the surface stress T can be brought arbitrarily close to zero. 

While cutting Kt, care is needed with surfaces where the bulk stress is singular. The neighbor- 
hood of a point on such surfaces is to be cut such that the length dimension parallel to the normal 
(of the surface) is arbitrarily small compare to other length dimensions. This way we will be left 
essentially with areas to be relaxed from stress, if any. The resulting stress-free configurations at 
the singular interface are of dimension one less than those obtained from the bulk. Their union 
forms a two dimensional smooth manifold N . A local configuration in J\f is identified with the 
local tangent space of N . If the tangent space T St (x) is mapped (locally) into two disjoint local 
configurations in J\f, for reasons that will become clear below, then we call the singular interface 
incoherent (at x). We denote the two local configurations by Tjj and Tjy- (in rest of the paper, a 
superscript 7 will represent an association with Tjj- configuration and 5 with T^; they are not to 
be confused as indices). Otherwise, if the mapping is injective then we call the singular interface 
coherent and denote the local configuration by Tfj. Incoherency of the interface implies a contin- 
uous distribution of dislocations over the interface; we postpone the discussion on this aspect till 
the next Section. The process of relaxation is illustrated through a cartoon in Figure 

Let K € Lin be the local map from tangent space Ki to V Kr at X € K r \ S r . Both H and 
K are assumed to be continuously differentiable except on the singular surface. The following 
decomposition 

H = FK VX G K r \ S r (2.78) 

is admitted (conventional plasticity theories usually represent tensors H and K _1 by F e and F p , 
respectively). Since we demand unloading to be elastic in nature, we call H the elastic distortion. 
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We identify K with plastic distortion, for reasons that will become apparent when we discuss 
dissipation in Subsection 14.11 Define distortion maps on the surface 

H 7 = H+l 7 , HP 5 = H~t 5 and (2.79) 
IK 7 = K + l 7 , K 5 = K~l 5 , (2.80) 



where superscripts 7 and S denote the association with the two local configurations in M at a fixed 

1%- and I s : V -»• T 5 N 

N 7 ® N 7 and I 5 = 1 - N* 5 ® N 5 , (2.81) 



material point. The projection tensors I 7 : V — >■ TX- and l* 5 : V — >• T£- are given by 



where N 7 and are unit normals to Tjj- and T^, respectively. They are related to N and 11 as (cf. 

m = |(H±)-^N«| ^ (2 - 82) 
™- K ' (2.83) 



|(K±)- T N«|' 

where a = {7, 5}. Here, superscript + appears with 7 and — with 5 (either pair can be equivalently 
used to obtain e or N). Normals N 7 and N 5 will coincide only if the jump [H _1 ]] (or [K _1 ]]) across 
the interface is of Hadamard's rank one form, i.e. if [H -1 ] = h(g>in (or [[K -1 ] = k(g)N for arbitrary 
h and k). Otherwise they will be distinct, resulting in two distinct local configurations after stress 
relaxation at each material point on the interface. The former case leads to a coherent interface 
and the latter to an incoherent interface. Let j a (jf) be the ratios of infinitesimal area on S r (st) 
to local infinitesimal areas on T^. They are given by 

f = KK^PH and f s = KH^N"! (2.84) 

and are related to each other as 

./(/') '• (2-85) 

where j is the ratio of infinitesimal area on st to S r . 

A straightforward calculation, using (|2.82p and (I2.83p . confirms that 

IF = lM a , and K a = tK a . (2.86) 
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(i) Incoherent interface without bulk dislocations 
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(ii) Incoherent interface with bulk dislocations. 
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Figure 3: A cartoon illustrating relaxed configurations. The rectangular blocks represent infinites- 
imal neighborhood in the body. In case (i) the bulk is free of defects and therefore the relaxed 
manifold consists of two disjoint global configurations which do not fit together at the interface. 
When the bulk is dislocated, as in case (ii), the relaxed manifold is a collection of only local 
configurations which do not fit together with each other. 
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Therefore HP map local configurations in N to T St and HP map local configurations in J\f to T$ r ■ 
Employing flU^D , (f!T80l) . ([2786]) . and (l23T]) i. and taking limiting values of (p778|) . we are led to the 
following multiplicative decompositions on the singular surface: (compare with (12.78P ) 

1° = ¥K a VX G S r , with a = {7, 5}. (2.87) 

There exist unique pseudoinverse tensors (HP) -1 and (HP) -1 such that (here and elsewhere, no 
summation is implied for repeated superscript a unless explicitly stated) 

(HP)" 1 HP = l a , HP (HI")- 1 = 1, (2.88) 

(IP)" 1 K a = l a , and K a (K Q )^ = 1. (2.89) 

Identities 

(HP) -1 = (H ± )~ 1 i and (MP) -1 = (K^) -1 1 (2.90) 

are then immediate, as can be verified by direct substitution using (|2.79p . (|2.80p . and (|2.86p . Here 
superscript + is used to define (HI 7 ) -1 (and (HP) -1 ) and — for (HP 5 ) 1 (and (H? 5 ) 

Finally, on the basis of our earlier discussion regarding F, we can obtain the following results 
for interfacial elastic and plastic distortions: (compare with (|2.34p . ()2.44p . and (|2.45p ) 

(IP)* = f (N ® N a ), j£ a = j a (K a )- T , = 0, j a = ri^K")- 1 • 1, (2.91) 
(M a Y = j s Q (n ® N a ), (jf) Ba = j«(M a y T , (jf) Na = 0, and jf = j^^y 1 ■ 1 (2.92) 

for each a = {7, 8}. 

2.5 Symmetries of space and matter 

To develop physically consistent constitutive models we need to exploit spatial and material sym- 
metries as afforded by the structure of the space and the material. We will obtain restrictions 
on the form of constitutive response functions (for example stress, energy, kinetic laws), both in 
bulk and on interface, upon imposing their invariance under the symmetries. In this subsection we 
discuss kinematical changes induced by the symmetries but postpone their application to response 
functions until Section [H 
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Material frame indifference The material response is expected to remain invariant under 
arbitrary changes in the frame of reference. For the case at hand this amounts to requiring invariance 
under superimposed rigid body rotations. The motion transforms as x — > Q(t)x + c(t), where 
Q € Orth + and c £ V are functions of time only (arrow denotes the transformation), and the 
deformation gradient as F — > QF. We assume T — > QTQ T for Cauchy stress and therefore 



The transformations for H and K, obtained by utilizing the constitutive framework of elastic 
plastic deformation in [16], are given by H —¥ QH and K —¥ K; while for unit normal in it follows 
from (I2.28p i . 11 —> Qe. All other unit normals: N and N a are invariant, where a = {"f,5}. As a 
consequence I ->■ QlQ T , 1^1, and t a -)• t a . Using these with (127791) and (127501) . we obtain 



HP -> QHP and K a -> K a . Moreover, (HP)" 1 -> (HP)" 1 Q T and -> (IK")" 1 , cf. (12790]) . 



Compatible changes in the reference configuration Our choice of reference configuration is 
arbitrary as long as it is a connected subset of the Euclidean space. We define compatible changes 
in the reference configuration as those maps which preserve its connectedness and Euclidean nature. 
A constitutive response function should be invariant with respect to arbitrary compatible changes 
in the reference configuration unless it is defined explicitly with respect to a particular reference 
configuration (which is the case in conventional elasticity theories, where the reference configuration 
is usually the stress- free configuration) [llj, \12\ HE] . 

At a fixed time consider a variation of K r from K ri to k T2 defined by a continuous but piecewise 
differentiable bijective map X2 = A(Xi), where Xi € K ri and X2 6 K r2 , with invertible gradient 
A = ViA, where Vi is the gradient with respect to Xi. To maintain compatibility of A at 
the singular surface, tensor A can at most have Hadamard's rank one jump discontinuity at the 
interface, i.e. [A] = d (g) Ni for arbitrary d € V. Local configurations in M. and Af, and the 
global spatial configuration ttt, all remain invariant under compatible changes in K r . Using obvious 
notation we have K^dXi = Kg 1 dX 2 away from the interface, leading to 



P -> QP. 



K 2 = AKi, H 2 = Hi, and F 2 = F]A 



(2.93) 



The unit normal Ni _L T$ ri at some Xi € S n changes to 




(2.94) 
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such that N2 _L Tg r at X2 = A(Xi) € S r2 , where ja is the ratio of infinitesimal areas on S ra and 
S ri . Other unit normal fields, N a (a = {7, 5}) and e remain invariant. For interfacial distortions, 
recall 1|277S]> . (jZBQD , and (^86j) to write 

Kg = AI" and Eg = Iff, (2.95) 

where A = A li. Moreover 

(K")^ 1 = (K )" 1 A -1 , (H")- 1 = (HP)^ 1 , and F 2 = FiA" 1 , (2.96) 

where A -1 is the pseudoinverse of A such that A~ 1 A = li, AA _1 = I2, and A" 1 = (A ± ) _1 I2, cf. 
([235]) and (|236|) . Relations (12^96]) can be verified by direct substitution using (pEggj) . ([2T89]) . and 

In fact, it turns out that A has to be continuous across the interface; i.e., d = 0. To show this 
we start by noting the change in U\, 

U 2 =Ja 1 Ja±U 1 . (2.97) 

This follows upon recalling the discussion leading to (|2.28j) and using | (A ± ) T Ni| = JaJa± from 
(I2.94p . For U2 to be single valued we require J4+ = J a-- Applying this to the identity Ja+ = 
J A -{\ + (A-^d • Ni), cf. Equation (1.1.6) in [38], we get d • (A")- T Ni = i.e., d G T Sr2 . A 
similar argument yields e £ Tg ri , where e satisfy [A -1 ] = e ® Next, consider an arbitrary 
single valued scalar field, ir, defined on S ri . Let ir be invariant under the transformation. Denote 
the normal time derivative of it with respect to S r2 by n. Use (|2.94p . (|2.97p . and the definition of 
normal time derivative to get 7r = tt + L^Vfvr • li(A =t )~ 1 N2. Subtract these two equations (one 
with superscript + and the other with superscript — ) from each other, and use the continuity of 7r, 
tt, and U2V1 tt, to deduce lie = or, since e € T$ ri , e = 0. Therefore, A -1 (or equivalently A) is 
continuous across the singular interface. 

We need to make two assumptions regarding the nature of A at the interface. We assume the 
normal time derivative to remain invariant i.e., tt = tt for arbitrary tt and A = 0. Both of these are 
motivated by the requirement of dissipation, given in (|4.18p below, to remain invariant under the 
transformation in reference configuration. The former of these assumptions leads to liA _1 N2 = 
implying 

A" 1 = /1N1 <g> N 2 + A" 1 (2.98) 
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at the interface, where h = JaJ^ ■ The second assumption, on using the definition of normal time 
derivative, yields 

V1AN1 = 0. (2.99) 
Finally, we evaluate the transformation of the curvature tensor L. Use (|2.98|) to rewrite (|2,94D 

as 

N 2 = frANi. (2.100) 
Differentiate it and apply (|2.99j) . after noting the symmetry in the gradient of A, to get 

L 2 = ZiALiA" 1 . (2.101) 

Material symmetry The concept of material symmetry is related to local configurations in the 
relaxed manifold A4. Let Q be the symmetry group associated with local configuration «j. An 
element G € Q then brings about a local change in Ki, modifying distortions from H and K to HG 
and KG, respectively, such that the constitutive response functions remain invariant. For solids, 
Q C Orth + . Consider, for example, the strain energy density W(H) (defined with respect to Ad). 
Then, W(H) = W(HG). We assume the material to be uniform away from the interface thereby 
allowing G to be discontinuous at the interface. Note that the material neighborhoods separated 
by an interface can have distinct symmetry groups, as is generally the case with poly-crystalline 
materials. 

The situation at the interface is more involved. For an incoherent interface we have two disjoint 
local configurations at each material point on S r and each of these can have a distinct symmetry 
group associated with it. Accordingly, we denote Q 1 and Q s as distinct symmetry groups associated 
with local configurations Tjj- and T^, respectively. Elements G 7 6 and G s € Q s thereby modify 
distortions H ± and K ± to H ± G Q and K ± G a , with a = 7 when the superscript is + and 5 when 
it is — (no summation implied), such that constitutive response functions on the interface remain 
invariant. Assume Q a C Orth + . This is sufficient to ensure that symmetry transformations neither 
strain the interface nor bring about a change in its local area, as shown below. Under the action 
of a symmetry map the interfacial normals N 7 and N 5 are modified to (G 7 ) T N 7 and (G 5 ) T N 5 , 
respectively. This follows immediately on using Nanson's formula and G a € Orth + . As a result 
identity tensors 1 Q transform to l a ' = (G a ) T l a G a , where l a ' represent surface identity tensors 
for transformed local configurations at the interface. Using these we can obtain the transformation 
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for interface distortion tensors, denned in (|2.79p and (|2.80p . as 



HP -> HPG Q and W 



•a 



a 



G' 



Q 



(2.102) 



where G 1 



Q 



l a G a (no summation over a). Tensors G a satisfy 



G a (G a ) T = l a , (G a ) T G a = l a ', and |(G a )*l = 1. 



(2.103) 



The last of these can be proved using (G a )* = G" and (1 Q )* = N a <8> N° (no summation on 
a). Relations (|2.103j) '? 3 ensure that symmetry maps do not introduce any additional strain at the 
interface and bring about no change in its local area. The former of these claims follows upon 
noting the definition of interfacial strains, 



where C a = (HP) HP. The latter claim is a consequence of the fact that |(G°)*| represents the 
change in the area of local configurations under the action of symmetry map. 

Remark 2.4. (Continuous and discrete symmetry groups) If the symmetry group Q is continuous 
(for example isotropy or transversely isotropy) then the dislocation density a (defined in f|3.22j) 
below) fails to be a characteristic of a body (cf. Theorem 8 in [36] and §1.2.5 in [I2])- Indeed, 
at a fixed time for a material point away from the interface, any two intermediate configurations 
and Ki 2 are related by K2(X) = Ki(X)G(X)A (cf. Theorem 2 in [36]), where G(X) £ Q 
and A € InvLin is constant. If the symmetry group is discrete then K2(X) can be a smooth 
field only if G is constant. However, for the continuous symmetry group, G can be a continuous 
function of X. This leads to a uniquely defined dislocation density tensor only in the former case 
(for details see the cited references). It should be noted that for metals such continuous symmetries 
are, in practice, used to model a poly-crystalline material which, at a sufficient macroscopic scale, 
is considered as a random aggregate of single crystals. To model a poly-crystalline material within 
our framework, we would need to construct the theory at the level of an individual grain boundary, 
which separates single crystals. This is precisely one of the motivations for the current work. 

3 Continuum distribution of dislocations: bulk and interface 

A distortion field (F, H, or K) is said to be compatible, away from the singular interface, if it is 
given by the gradient of a differentiable vector field; and compatible at the interface if it admits a 



a 




2 



(2.104) 
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rank one jump of Hadamard's type across the interface. Within our model we assume both H and 
K to be incompatible, and F to be compatible, everywhere in the body at all times. As we shall 
see below in Subsection 13. 1\ incompatibility of H (or equivalently of K) in effect contributes to 
the Burgers vector for arbitrary circuits in m (or K r ) thus affording a relationship with continuous 
distributions of dislocations. The vanishing of Burgers vector is equivalent to the compatibility of 
distortion fields H and K in the region. However, caution is required when interpreting a continuous 
distribution of dislocations in terms of a distribution of discrete dislocations. It is possible to have 
multiple arrangements of discrete dislocations which correspond to the same continuous dislocation 
density. This point, in the context of interface dislocation density, has been well illustrated in 
Section 2.4 of [43]. Moreover, interface dislocation density can be represented by a two-dimensional 
array of discrete dislocations only for low magnitudes of dislocation density. This follows from the 
fact that individual dislocations cannot be placed arbitrarily close to each other (cf. Section 38 in 
[10j). Similar remarks hold for bulk dislocation density. 

The dislocation densities obtained in Subsection 13. II are not invariant under compatible changes 
in the reference configuration. This is unacceptable for a genuine measure of defect content which 
should be independent of the choice of a reference configuration. We are led to obtain invariant 
measures of bulk and interface dislocation densities which thus qualify as argument of constitutive 
response functions associated with plastic evolution, cf. [12], [16] for bulk and Subsections 14.11 and 
14.21 for interface. 

The dislocation densities associated with bulk and interface are not independent of each other. 
As shown in Subsection 13.31 the projection of the jump in bulk dislocation density along the 
interface normal is related to the (surface) divergence of interface dislocation density. This relation 
immediately proves that interface dislocation density, in contrast to its bulk counterpart, does not 
have a vanishing divergence. 

3.1 Interface dislocation density 

Define Burgers vector 



where C is a close material curve which intersects S r in finite number of points and K is piecewise 
smooth. The plastic distortion is singular on S r and therefore the integral in (13. ID will have 




(3.1) 
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singularities only over a set of measure zero (a finite collection of points on a line constitute such a 
set). Let Aq be the area of the surface enclosed by C and let T = Aq H S r be the curve which lies 
at the intersection of Ac and S r . Use the Stokes' theorem for piecewise continuous tensor fields 
(see for example [16]) to obtain 

I K" 1 dX= f (CurlKrYNcdA- [ {K^jdX, (3.2) 
J c J a c Jr 

where Nc* is the unit normal field associated with Ac- Let {ti, t2} € Ts r (X) be such that {ti, t2, N} 

is a positively oriented orthogonal basis at X £ S r . Orient t2 along the curve T so that dX = 

t2<iL = (N x ti)dL on T. Define a tensor field (3 r on 5 r , the (referential) interface dislocation 

density, such that 

[K^Kti x N) = $rt\. (3.3) 

Recognizing CurlK -1 as the (referential) bulk dislocation density, denoted by a r , the net Burgers 
vector associated with C can be written as a function of dislocation densities: 

B(C,t)= / aJ~N c dA+ [ $ T r t x dL. (3.4) 
J a c Jr 

Writing [[K" 1 ]] = [[K _1 ]]l with 1 = N ® N + 1 ® t a (a = 1, 2), §2M leads to 

[K- 1 ] =k®N-/fe (N) , (3.5) 

where k € V is arbitrary and 

e (N) = *1 ® *2 - *2 <8> ti (3.6) 

is the two dimensional permutation tensor density on Ts r (X). It satisfies emj = Re^)R T for all 
two dimensional orthogonal transformations R that preserve the orientation of Ts r at a material 
point. Therefore, any pair of vectors in Ts r (X) which with N form a positively oriented orthogonal 
basis may be used in the definition of e^y 

Projecting (|3.5p onto 1 and applying (|2.90p 9 yields 

[K- 1 ] 1 = (IPT 1 - (K 5 )- 1 = -fis in) (3.7) 

which, on using = — 1, leads to 

0?1 = [K- 1 ]^. (3.8) 
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This determines the action of f3j on Ts r (X). The action of 0^ on N is indeterminate and may be 
set to zero without loss of generality. Consequently we assume 0£ to be superficial (i.e., 0^1 = f3 r 



wiuiuuL iuss ui geneiciiiiy. ^unsequeniiy we assume 
and write 

>9r = PK _1 I^CN)- ( 3 - 9 ) 
The interface S r is coherent at X € S r if the Burgers vector as defined in (13. 4h has no contri- 
bution from the line integral for all closed curves C such that X G T. Accordingly, SV is coherent 
at X € <SV if and only if the interface dislocation density j3 r at X vanishes. 

Proceeding in parallel we can derive the spatial form of above relations. In particular 

[H- 1 ] = h ® e - 0f e (n) and = p" 1 ^^), (3.10) 

where h G V is arbitrary and is the two dimensional permutation tensor density on T St (x). 
The tensor /3 t is the (spatial) interface dislocation density which satisfies 0fl = 0f . The Burgers 
vector using spatial description is given by 

b(c,t) = / H _1 dx = / (curlH~ 1 ) T n c da - /pr^dx 

= / (a t fn c da+ ! 0£i 1 dl, (3.11) 
J A c J-y 

where c is a (closed) spatial curve enclosing area A c (with unit normal n c ) and 7 = A c n St- The 
tensor at is to be identified with the (spatial) bulk dislocation density. The curve 7 is parametrized 
by arc-length I and has an associated tangent vector £2 such that the triad {ti,t2, m} forms a 
positively oriented orthonormal basis on St- 

The referential and spatial interface dislocation densities are not independent. To obtain the 
relation start by noting that the jump in F _1 across the singular surfaces has the Hadamard's form 
[[F" 1 ] = a ® e with a G V arbitrary. Using this and (|2.78p together with 

[H- 1 ! = (K- 1 )^" 1 ! + [K- 1 ]^- 1 ) (3.12) 

we derive 

h <g> e - 01 e (n) = (K- X )a <g> e + k <g> (F~ r )N - (3js {N) (F^ 1 ). (3.13) 

Nanson's formula (|2.28P i ensures that (F~ T )N is parallel to n. Multiplying (|3. 13|) on the right by 
E( n \ then furnishes 0[ in terms of f3 r , 

(3? = -/3^ (N) (F- 1 )e (n)! (3.14) 
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while the normal component of (|3.13p yields a relationship among a, k, and h: 

h = (K-^a + (m • (F~ r )N)k - 0?e w (F _1 )ii. (3.15) 

Remark 3.1. If K~ x is the gradient of a piecewise continuously differentiable deformation i.e. 
K 1 = Vx p , then 05] 

[K^ 1 ] = k (g) N + V s |x p l- (3-16) 
The interface dislocation density, given by (|3.8p . then has the form 

fit = VVle(N) = V[f]e W , (3.17) 

where we have used the identity lem) = £(n)- This situation occurs when the bulk dislocation 
density at r = CurlK" 1 vanishes on either side of the interface and the dislocation distribution is 
restricted to the singular surface S r . An equivalent formulation holds with respect to the spatial 
configuration. The previous work on incoherent interfaces [6j [33] has been in fact restricted to this 
case with the exception of the paper by Cermelli and Sellers [9] who have developed the theory in 
the context of crystal lattice vectors. 

Remark 3.2. (Interface dislocation nodes) An interface dislocation node, as introduced by Bilby 
[2], is the line of intersection of interfaces with dislocation density distributions. The analysis in 
[2] is restricted to plane interfaces and infinitesimal strains. We extend it for curved surfaces and 
finite distortions. Consider iV surfaces intersecting at a line L C K r . Each surface has an associated 
normal and a distribution of interface dislocation density. The following compatibility relation 
holds in a neighborhood infinitesimal close to L 

N 

^pC" 1 ®] = 0, (3.18) 

i=l 

where the index i in the superscript represents the i'th interface. This relation follows from the 
observation that on passing around the line L (in a small neighborhood) one reaches the initial 
material point; and (K -1 W) + = (K^ 1 ^ 1 )) - , where (i + 1) should be taken as 1 when i = N. 
Use K33J) to rewrite (pUBD as 

N 

^(k«®N»-/3j% (l)) ) = 0. (3.19) 
i=i 
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Let t be the unit tangent vector field associated with line L. We can therefore choose vector 
€ Tg(i) at X € L such that {t, t^,N^'} forms a positively oriented orthonormal basis at X G L 
for each intersecting surface. We also have, cf. (j3.6|) . 

£(NW) = * ® t (<) - ® t. (3.20) 
On substituting this in f|3. 19j) it follows immediately that 

N N 

= ^[K-^Jt = ^/3^t« (3.21) 

i=l i=l 

at X € L. The outer equality in (|3.2ip provides us with a compatibility condition relating the 
interface dislocation density tensors of various intersecting surfaces. This can be compared to the 

equation of conservation of Burgers vectors. Equivalently, in terms of the spatial surface dislocation 

N TO- 

density, we can obtain Y,Pt * W = °> where * W e T (*)■ 

i=i St 

3.2 True interface dislocation density 

A measure of the bulk dislocation density, invariant with respect to compatible changes in the 
reference configuration, is given by [8] 

JxKT 1 CurlKT 1 = a = J^H" 1 curlH" 1 . (3.22) 

Tensor a. is a map from a local configuration in M. \ M onto itself. It is, therefore, also invariant 
under superimposed rigid body motions. The importance of a is perhaps most evident in its 
appearance in constitutive response functions related to plastic flow (cf. [HI [16] and (|4.68p below). 
Response functions, if invariant under compatible changes in the reference configuration, cannot 
depend explicitly on K and their dependence on VK is only through a [HI [TT] . 

We now obtain an analogous measure for interface dislocation density. Consider two reference 
configurations related by a compatible deformation with smooth gradient A (refer to Subsection [23] 
for the kinematics of a compatible change in the reference configuration). We impose the require- 
ment that under a compatible transformation the Burgers vector is left invariant and consequently 

/ (Curl! K^fVfadA!- [ pC^JMXi = / (Curl 2 K^) T N C2 dA 2 - [ lK^}dX 2 , (3.23) 
Ja Ci JVi JAc 2 Jt 2 



32 



where A C2 = X(A Cl ), T 2 = A(ri), and K 2 = AKi (cf. fZMb ). The area integrals in (|3^3|) are 
equal, owing to the invariance of a, thereby reducing (|3.23p to 

/ [Kr^dXx = [ iK^jdXz. (3.24) 

Let t a and L a be the unit tangent vector and the arc- length associated with r a , respectively 
(a = 1,2). Then 

t 2 dL 2 = dX 2 = AdXi = AtidLi. (3.25) 
Substitute dX. 2 from ()3.25|) into (|3.24|) and employ the arbitrariness of T\ to obtain 

[K^JAti = [Kf 1 ]*! (3.26) 

for ti € Ts ri ■ Because tii is otherwise arbitrary this relation is satisfied only if 

[K^A-lK^Hc^Ni, (3.27) 

where c € V is arbitrary. Project both sides of (|3.2T|> on li and use (|3.7p to get 

tf 2 e m A = (3^s iNl) (3.28) 

or equivalently, on applying (|2.95|) i . 

^e (Na) K§' = ^ 1 e (Nl) K? (3.29) 

with a = {7,5}. This leads us to define invariant (or true) interface dislocation densities /3 7 and 
(3 S by 

pn = ^£ (N) K 7 and (3 5 = ff £{N) K 5 . (3.30) 
Define the relative distortion tensor 

M = (H 7 )~ 1 IHI <5 = (IC 7 ) _1 1K <5 (3.31) 

which is a linear map between two local configurations at a material point on the incoherent 
interface. It was introduced by Ceremelli and Gurtin [61 [7] where it was called the incoherency 
tensor. The second equality in (|3,31|) follows from (|2.87|) . It is checked easily that M satisfies 
1 5 and M = 1 7 M. Also, observe that 

/3 5 = /3 7 M. (3.32) 
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Substitute 0£e(K\ from (|3.7p 9 and use (|3.3ip to get 



/3 7 = M 



I 7 and (3 s = I s - M, 



(3.33) 



where M" 1 is the pseudoinverse of M such that M _1 M = I s and MM -1 = I 7 . 

If in the above analysis we substitute F for A we obtain, instead of (|3.29p . the following relation 



with a = {7, 5}. The true interface dislocation densities can thus be equivalently expressed in 
terms of elastic distortion. 

For the incoherent interface the tangent plane (to the singular surface) in the reference (or 
spatial) configuration is mapped (locally) into two tangent planes in the relaxed manifold. As a 
result we have two measures, /3 7 and /3 5 , of true interface dislocation density for each X € T$ r (or 
x £ T St ). This is in contrast to the bulk where, as pointed out in the beginning of this Subsection, 
we have a single measure of invariant bulk dislocation density. 

Remark 3.3. (Equivalence between different dislocation densities) With different measures of dislo- 
cation density distributions, both in bulk and on interface, it is useful to investigate the possibility 
of equivalence between them. The true bulk dislocation density, a, vanishes if and only if a r 
(or at) vanishes. This is evident from (|3.22p . where the determinants are positive and distortions 
invertible. This however is not the case with their time derivatives. Similarly, the true interface 
dislocation densities vanish if and only if (3 r (or f3 t ) vanishes. This follows on using (I2.88P and 
(|2.89p in (|3.30p and (|3.34[) . Their normal time derivative however might not all be zero at the same 
instant. 

3.3 Relationship between bulk and interface dislocation densities 

We now relate the jump in bulk dislocation density across the interface with interface dislocation 
density. An integral form of conservation of dislocations, for an arbitrary volume which intersects 
the singular interface, is also obtained. Our discussion involves referential dislocation densities; 
similar relations can be obtained using their spatial counterparts. 

A compatibility condition for VK" 1 , discontinuous across S r but otherwise smooth, is given by 



Pf e {n) U a = Pie w K 



a 



(3.34) 



!H|45] 



(3.35) 
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where Q € Lin is arbitrary. In terms of indicial notation this is alternatively written as 

[if-Jj = Q jt N k + iKr^l^. (3.36) 

Multiply (|3.36p throughout by e^i and use the definition of referential bulk dislocation density to 
obtain (the subscript present, if any, in the bold notation is written as a superscript in the indicial 
notation) 

Mjl = emQjiNk + emlK^j^tmk. (3.37) 
Using e^NiNfc = 0, the normal projection of (I3.37P yields 

KP.^.fcil^ 1 ],*. (3.38) 
The jump in K _1 , given in ()3.5|) . can be written in indicial notation as 

Replace {Kf. 1 } from (13391) in (ETggl) to get 

[o^-jNi = eifc,AyNj >fc Ni + e lkl (^ gj e^ ) ) :k N i . (3.40) 

The first term on the right hand side vanishes since ej/yN/^Nj = 0, which can be proved using (|2.6I) i 
and the skew symmetry of en-i- Consequently (j3.4Q[) reduces to 

[a<,]N< = ^.^weJ'Ni + ^ e «e£N, (3.41) 

We note the following two identities: 

e ik l^i = l qk and (3.42) 

eawe^Ni = 2FN 9 . (3.43) 

Relation (|3.42p follows from the definition of e^. The proof for (|3.43p . which is left to the reader, 
is however more involved and requires calculating divergence of (I3.42P and using VN = — L + 
(VN)N (g> N, where L is the symmetric curvature tensor defined in (|2.16p . 

Use ([5T12]) and (f3^3l) to write (l3T3Tj) equivalently as [o^JN = Div s / 3^ + 2H(3^N or 

[c$pi = Div s $ VX e S r , (3.44) 
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given that /3 r is superficial. This is an important result highlighting the nature of interface dislo- 
cation densities (cf. Div s a^T = VX G K r \S r ). To expand on this we use the surface divergence 
theorem (|2.47p . For an arbitrary surface 5 C S r we have 

/ 0^.udL = [ Biv s (3jdA, (3.45) 
JdS Js 

where v is the outer unit normal to dS. In addition, use the divergence theorem for piecewise 
smooth fields, and the identity Div = away from S, to derive 

/ exjNdA = / [aJT]NdA, (3.46) 
Jan Js 

where N is the normal to dCl. Combining (|3.44p . (|3.45p . and (|3.46|) yields 

/ fi T r vdL = [ aJlSidA (3.47) 
JdS Jan 

as the integral law for conservation of dislocations in an arbitrary volume ft C K r such that S = 
S r n fi / 0. 

If [oJTJN = 0, i.e. there are no external sources to interfacial dislocation density, then (|3.47p 
reduces to 

ffivdL = 0, (3.48) 



'as 

which can be interpreted as the conservation law for interface dislocations. This is analogous to 
the conservation law for a r according to which, for an arbitrary volume O C K r with S r H = 0, 

cx^NdA = 0. (3.49) 

an 

Relation (|3.49p imposes the restriction on bulk dislocations to not end arbitrarily inside CI. A par- 
allel interpretation in the context of interface dislocation densities is furnished by (|3.48p . Therefore 
for a vanishing normal jump in a£ an incoherent interface cannot end arbitrarily inside the solid. 
Such an interface will either end at the boundary of the solid or at a surface dislocation node (see 
Remark 13.21 above) . 



Remark 3.4. For moving interfaces (U ^ 0) we can derive an alternate expression for the jump in 
bulk dislocation density. The jump in KJ, 1 can be expressed in terms of the normal time derivative 
of \KJ?l as 

IK' 1 f= UlK-^Nk + \K-% (3.50) 
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which is obtained by subtracting the normal time derivative of yK-^J from that of yK-j 1 
Substituting (EOHj) into (EUKl yields 

UQji = 1^1-1X^1 (3.51) 

Replace Q in (I3.37D from ()3.5ip to obtain 

£%y = -eMlKr^k + ewlKJ? fl°N fc + U e^RT^l^ (3.52) 

as the jump condition for bulk dislocation density across a moving interface. This relation reveals 
various sources in the production of bulk dislocation density (near the interface) as the interface 
moves with velocity U . Observe that [o; r ] does not necessarily vanish for continuous K" 1 . 

4 Interfacial plasticity: dissipation, energetics, and kinetics 

Plastic flow is a dissipative process involving irreversible restructuring of the microstructure which 
in turn affects the macroscopic behavior of bodies. Away from the interface the dissipation is 
caused by evolution of plastic distortion. At the interface there are three dissipative mechanisms: 
the motion of interface (governed by its normal velocity), the evolution of plastic distortion, and 
the evolution of relative plastic distortion (incoherency) . All these will, in general, be coupled 
to each other. In the following we start with specific constitutive assumptions on the nature of 
energy densities and stresses and then use them to evaluate the dissipation. The driving forces for 
various dissipative mechanisms are obtained. Based on the list of dissipative fluxes and driving 
forces general forms for kinetic laws are proposed and then simplified using invariance requirements 
under various symmetries. 

4.1 Dissipation inequality 

We now revisit dissipation inequalities (|2.64p and (J2.70P assuming energy densities such that [TB] 

^ = J-V(H) and (4.1) 

$ = (j7)-i w (eT ) e <5 ), (4.2) 

where W is the bulk energy density per unit volume of a relaxed configuration and w is the 
interface energy per unit area of a local configuration Tjj- (we can equivalently consider an interfacial 
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energy per unit area of T^). While it is possible to include higher-order gradients at the interface 
to reflect bending and other weakly non-local effects (as in [42] ) . it is our view that the present 
'membrane- like' model, in the spirit of the Gurtin-Murdoch model of interfaces [21], represents the 
leading-order effects faithfully if the actual interfacial region is sufficiently thin (as distinct from 
our present representation as a discontinuity surface). This issue represents a direction for future 
research in interfacial plasticity. 

Under the hypothesis of hyperelastic response (during elastic unloading), the bulk and interfacial 
Cauchy stresses are assumed to be 

J H T = W H H- T and (4.3) 
j]T= ^(M a f, (4.4) 

a=7,5 

respectively, where the summation is over both the local configurations in A/" at a fixed material 
point on the interface. Recalling P = TF*, in addition to (j2"760l) . ([2785]) . and ([2787]) . leads to the 
corresponding Piola stresses 

J K P = W H K T and (4.5) 
fF = £ mm* (K a f • (4-6) 

a=7,5 

We note that relations (|4.4p and (|4.6|) are motivated by the assumption that the total interfacial 
Piola stress P is power- conjugate to interfacial deformation gradient F, as in (|4.13|) below. 

The dissipation inequality in the bulk (|2.64p . on using (|4.ip and (I4.3p . reduces to (see for 
example [T3] [T6]) 

E • KK 1 > VX G K r \ S r , (4.7) 

where E is the bulk Eshelby tensor defined in (12.7ip . 

We now evaluate interface dissipation inequality ()2.70p under the above constitutive assump- 
tions. Use (14. 2p to obtain the normal time derivative of interface energy density 

$ = -p a 7 ) -1 * + on -1 E ■ ( 4 - 8 ) 

o=7,5 

Additionally, note that i«iaN ft = (no summation) for a = {7, 5} (cf. (|2.42|) i ). Substituting p 
from ([2311 ? into (JO) then yields 

$ = — K 7 (IC 7 )~ 1 • 1$ + {py 1 w ma l a ■ BP. (4.9) 

a=7,<5 
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Taking the normal time derivative of (|2.87p and using it to replace EI Q above leads to 

$ = — K^K 7 )"" 1 • 1$ + (j 7 ) -1 ^2 (w ma (K a ) T ■¥ + ¥ T w ma t a -K a Y (4.10) 

a=7,<5 

or equivalently (recall (|4.6|) and (|2.89]) i) 

4» = — K 7 (K 7 ) _1 • 1$ + P • F + (j 7 )" 1 ^ (F T w ma (K a ) T ■ K a (K a )'' 1 ^ . (4.11) 

o=7,5 

This can be expressed succinctly on introducing 

j 7 P 7 = w m {K~i) T and j 7 P 5 = w M 6{K 5 ) T (4.12) 

as two components of the interface Piola stress such that P = P 7 + P" 5 . Equation (|4.11|) can be then 
rewritten as 

<|> = P . F — E 7 • i 7 ^ 7 ) -1 + ¥ T F S • K a (K 5 ) -1 , (4.13) 

where 

E 7 = $1-F T P 7 (4.14) 
is the interface Eshelby tensor associated with Tjj-. Define K s such that 

E = E 7 + E 5 , (4.15) 

where E is the total interface Eshelby tensor introduced in (|2.72p . Hence 



The apparent asymmetry in the definition of two interface Eshelby tensors, in (|4,14p and (|4.16p . 

r 7 



is due to our use of an interface energy density measured per unit area of Tlr. The dissipation 



inequality (I2.70p . on replacing $ from (I4.13p . acquires the form 

U ( N • [EJN + ^C/ 2 pr[|FN| 2 ] + E • L ) + ^ E Q • K^r)" 1 > VX e S r , (4.17) 



o=7,(5 

or 



where 



Ui + E ° ' ^(^y 1 > VX e S r , (4.18) 

a=7,<5 



N- [Ep+^C/V[|FN| 2 ] ) 
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is the driving force associated with the normal motion of the interface and L is the curvature 
tensor. The left side of inequality (|4.18p is the net dissipation caused by a moving interface and 
interfacial plastic flow. For a coherent interface the two local configurations Tjj and Tjy- coincide, 
and IK 7 = K s (= K, say). The dissipation inequality ([4.18P then takes the form 

[/f + E.K(I)^>0VXe5 r , (4.20) 

The interface Eshelby tensor appears naturally as the driving force for plastic flow at the 
interface. This is comparable to the plastic behavior away from the interface, cf. (14. 7p . where the 
bulk Eshelby tensor drives the bulk plastic flow. Moreover, it is clear from (|4.18p that the normal 
projections of plastic distortion rates, i.e. K a N a , do not participate in dissipation. This is not 
the case when the energy density depends explicitly on interface normals, as discussed briefly in 
Remark 14.21 below. 

Remark 4.1. (Area and volume preserving plastic flow) Assume Jk = 1 and j a = 1 for a = {7, 5}. 
Therefore, the plastic distortion brings about no change in the volume and the area of bulk and 
interface, respectively. The dissipation inequalities (|4.7p and (|4.18p are reduced to 

F T P • KKT 1 < VX G n r \ S r and (4.21) 
US- ^2 FTp ° ' > VX € S r . (4.22) 

a=7,5 

This form is similar to classical plasticity theories wherein the plastic flow is driven by stress rather 
than the Eshelby tensor. 

Remark 4.2. Consider an interfacial energy density of the form w(W, U s , N 7 , N 5 ) with explicit 
dependence on the interfacial normals. Such energies (for coherent and unstrained interface) have 
been used for example to model the anisotropy of the surface during crystal growth |24j . The 
normal projections ^H a N Q (for each a = {7, 5}) are no longer zero, cf. (|2.42|) i . and they contribute 
to the driving force for the normal evolution of plastic distortion. Indeed, the term 

(w M *n a ■ M a N a + w m ■ N a ) (4.23) 

a=7,<5 

has to be now appended to the left hand side of the inequality (|4.18p . Take normal time derivatives 
of H Q N Q = and K a N a = to obtain EPN a = -HPN a and N a = -(K )- 1 !^", respectively. 
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Using these we can rewrite f|4.23j) as 

€ Q • (K a ) _1 i a N a , (4.24) 

a=7,5 

where «P = (M a ) T WM a ^ a — wif, cf. paragraph before (|2.73p . 

Interface energy density Invariance of w under a superimposed rigid body rotation requires 

w (BP,m s ) = w(QHF, QM 5 ) (4.25) 

for arbitrary Q € Orth + (cf. Subsection I2.5p . Tensors HP (a = {7, 5}) admit polar decomposition 
HP = R a lLP, where R a € Orth are non- unique and positive semidefinite tensors CP £ Sym are 
unique (cf. flggfl and (f2~38l) 'k U a satisfy U° = l a IPl Q . Define R a = M Q (U a ) _1 , where (IP) -1 
is the pseudoinverse of HP such that (ILP)" 1 HP = HP (HP)" 1 = 1 Q . Consequently relations IP = 
R a t a and QR Q ) -1 = (M a ) T hold, where (W*)' 1 is the pseudoinverse of R a satisfying (M a ) _1 ]R a = 
1 Q and R a (M a ) _1 = 1. Therefore, we also have BP (KP) T = I and (KP) T ir = t a . To this end, 
note the decomposition of interface elastic distortion tensors 

HP = 1PU Q (4.26) 

into two unique tensors. 

Choose Q = (W) T + N 7 <g> n; thus det Q = detR 7 (R 7 N 7 • 11) = (detR 7 ) 2 (= 1), where the first 
equality follows upon substituting W = R 7 — R 7 N 7 <8> N 7 in the expression for Q and then using 
identity (1.1.6) from [38] . To prove the second equality, recall (|2.92p i to write 

M 7 U 7 ti (g) M 7 U 7 t 2 = |(H 7 )*N 7 |m (4.27) 

and let t a € Tjj- (a = {1, 2}) be the two principal vectors of U 7 (the third one is given by N 7 ) such 
that {ti,t2,N 7 } forms a positively oriented orthogonal basis. The details are left to the reader. 
Substitute the assumed Q in (|4.25p to obtain 

w = w{W, (W) T M. 5 ) (4.28) 

or equivalently 

w = tD(C 7 ,M), (4.29) 
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where C 7 = (U 7 ) 2 and M = (HI 7 )" 1 !! 5 (= (IK 7 )" 1 !!' 5 ) is the relative elastic (or plastic) distortion 
between two local configurations at a material point on the incoherent interface, cf. (|3.3ip . Tensor 
M is a linear map from to Tfy and is related to true interface dislocation densities, as shown in 

Additional insight is furnished by rewriting the dissipation inequality (|4.18|) with the interfacial 
energy density given by (14. 29ft . Recall (14. 4p and define, cf. (14. 12f) . 

j 7 T 7 = w m {W) T and j]T 5 = w m s {M 5 ) T (4.30) 

such that T = T 7 + T s . The Cauchy stress tensor T is symmetric (cf. ()2.62p 9) and satisfies 
T = 1T1. Unlike T, T a (a = {7, 5}) are not symmetric. Their asymmetric parts are related as 
Skw(T^) = —Skw(T s ) (since Skw(T) = 0), where Skw() denotes the skew-symmetric part of the 
tensor (similarly let Sym() denote the symmetric part of the tensor). Noting that wm = 1 7 ^m1 <5 
and w c ~f = l 7 ?^! 7 (cf. (g32D i) we have, from K29\\ and (pT31~)) i . 

w = Sym(wc~i) ■ C 7 + % • M 

= (2WSym(iBci) - (W)- T w m M T ) ■ W + {W)~ T w m ■ M s . (4.31) 

The coefficient of HP above should be equal to ww* for each a = {7, 5}. Exploiting this correspon- 
dence and using (|4.30p we get 

j]T = 2IFSym(wcj)(Mr') T and j]T s = (M 7 )- T u) M (]H <5 ) T . (4.32) 

Therefore, the dependence of w on C 7 alone contributes to the total interface stress; and the 
dependence on M to the dissipation, as shown below. Use (I2.60p . (I4.6p . and ()4.12p 9 to obtain 

fW = 2WSym(w c ~ / )0&t) T and f¥ s = (EPy T w M 0& S ) T (4.33) 

for interface Piola stresses. On the other hand, use the normal time derivative of ()3.3ip 9 and the 
definition (|4.16p for E 5 to show 

E 5 ■ K^K 6 )" 1 = E 5 ■ (K^MiK 5 )- 1 + IK 7 ^ 7 )" 1 ) (4.34) 

and thus 

E Q • K ^ )- 1 = E • IK 7 ^ 7 )" 1 - (i 7 )" 1 ^ • M, (4.35) 

a=7,<5 
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where we have used (|4.16p . (|4.33|) -?. and (|2.87p . Substitute (|4.35p into the dissipation inequality 
(|4.18p to obtain its alternate form 

Uf + E • ^(K 7 ) -1 - if)- 1 tins. ■ M > VX € S r , (4.36) 

where f and E are given in (I4.19P and (I2.72P , respectively. Hence there are three dissipative processes 
active at an incoherent interface: the normal motion of the interface (driven by f), the evolution of 
plastic distortion (driven by E), and the evolution of relative plastic (or elastic) distortion (driven 
by (j 7 ) _1- u)m). The normal projections ]K 7 N 7 and MN 5 do not contribute to the dissipation. If, 
however, the energy density depends explicitly on interfacial normals then these normal projections 
would also contribute to the dissipation, as discussed in a remark above. The dissipation inequality 
as given in (|4.18p or (|4.36|) has to be satisfied by all kinetic laws governing the irreversible process 
of coupled plastic flow for an interface moving within a solid. 

Next, we note the restrictions imposed by material symmetry on the form of interface energy 
density. Recall the discussion in Subsection 12.51 and consider G 7 € £/ 7 and G 5 £ Q s , where £/ 7 and 

are the symmetry groups associated with Tjj- and T^, respectively. The material response at 
the interface remains invariant under the action of these groups. Therefore 

w{W,U 5 ) = w(H 7 G 7 ,IEff 5 G' 5 ), (4.37) 

where G a = t a G a , cf. (|2.102|) and (|2. 1U3|) . For the energy density given by (|4.29|) we require 

tZ)(C 7 , M) = w((G 7 ) T C 7 G 7 , (G 7 ) T MG 5 ) (4.38) 

to be satisfied by all elements of symmetry groups and Q 5 . Indeed, under the action of these 
symmetry groups C 7 = (M 7 ) T EI 7 transforms to (G 7 ) T C 7 G 7 , cf. ([2TT02]k and M = (M 7 )~ 1 IHI <5 trans- 
forms to (G 7 ) _1 MG 5 , where (G 7 ) -1 is the pseudoinverse of G 7 satisfying (G 7 ) -1 G 7 = I 7 and 
G 7 (G 7 )" 1 = I 7 (here I 7 ' represents the surface identity tensor for the transformed local config- 
uration at the 7-surface). Finally, recall (|2.103p and the uniqueness of the pseudoinverse to note 
that (G 7 ) -1 = (G 7 ) T . These considerations provide a point of departure for the phenomenologi- 
cal modeling of the behavior of incoherent interfaces which, contrary to coherent boundaries, are 
affected by two possibly distinct symmetry groups. 

Remark 4.3. The skew components of T a are conjugate to the relative spin tensor across the 
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interface. To elaborate, we write the normal time derivative of w in the form (using ([4 .30 ft ) 

j]w= ja - L ^ ( 4 - 39 ) 

a=7,<5 

where L a = H a (E[ a ) _1 . Decompose L a into symmetric (D Q ) and skew (W a ) parts. The tensors 
D a and W a are identified as the stretching tensor and the spin tensor, respectively. Substitute the 
decomposition in (I4.39h to get 

j]w = Yl (Sym(T a ) • D° + Skw(T a ) • W Q ) 

a=7,<5 

= Sym{T a ) • D Q + Skw(T 6 ) ■ (w 5 - W 7 ) , (4.40) 

ct=7,<5 

where the second equality is a result of the symmetry of T. Finally, note that Skw(T' y ) is completely 
determined from Skvu(wya), cf. (I4.32|) 9. 

Remark 4.4. If the interface is unstrained, i.e. HP = M a (for {a = 7, 6}), then the representation 
(p~29|) reduces to 

w = u)(l 7 ,M), (4.41) 

where K = (M 7 )' r ]R <5 is the relative rotation at the interface. Such energies are widely studied in 
the literature on metal interfaces, in particular grain boundaries, and many experimental methods 
have been devised for their evaluation. 

4.2 Kinetic relations 

Motivated by the dissipation inequality (|4.18p we consider constitutive functions of the form 

/ = /(ff,H^IK a ,M Q ,i Q ,]L,N a ,?7), (4.42) 

where the arguments of / can include both 7 and 5 variables (a = {7, <5}). Relations of the type 
/ = would then furnish possible candidates for kinetic laws at the interface. These relations 
complete the set of equations necessary to determine the evolution of state variables. 

We first obtain restrictions on / to be invariant under compatible changes in the reference 
configuration. Recall that a compatible change ensures that the Euclidean nature of the reference 
configuration remains unaltered. Consider, as in Subsections 12.51 and 13.21 two reference configura- 
tions K n and K r2 related by a map A such that X2 = A(Xi), where Xi £ K n and X2 G K r2 , with 
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continuous A £ InvLin given by A = ViA. As before, we assume that the transformation leaves 
the normal time derivative invariant and that A = 0. They imply conditions ()2.98j) and (|2.99p . 
respectively. 

The function / is invariant under a change of reference configuration from 

/(ffi, Uf, Kf , Hf , K?,Li, Nf , C/i) = /(ff 2> H£, K? , , K^,L 2) Nf , c/ 2 ), (4.43) 

where all normal time derivatives are given with respect to ft n . 

The transformations for bulk and interface distortions are given in (|2.93|) and (|2.95|) ; and for 
the normal velocity, normal, and curvature tensor in (|2.97p . (|2.100p . and (|2.10ip . respectively. 
Substitute A = K^(K^) -1 and ja = i 2 (if) 1 ( no summation over a = {7,/3}) into (|2.97p and 
(IXTUTp to get J2«/~±^2 = 3i J x± u i and i2 _lj it±( K 2) _lL 2 K 2 = if _1 ^±( K i) -1 LiKf; thus define 

U a = h a U and L a = (fc a )- 1 (K a ) _1 LX a , (4.44) 

where /i a = j a J~l_ (superscript + is used for defining /i 7 and — for h & ), as invariant (or true) 
normal speeds and invariant (or true) curvature tensors associated with the singular interface. 

The transformation in the driving force f can be evaluated by recalling its expression from 
(|4.19p and noting that 

E 2 = J~ A 1 A" T Ei A T , Pr2 = JjVri, and (4.45) 
E 2 = j~ 1 A- T EiA T . (4.46) 

These can be proved using OTTTD . (I2T72D . (jlTB) . (j43|) . and (p^3l) - (p^6|) with the fact that * and 
p r are densities per unit volume of K r and $ is a density per unit area of S r . Combine them with 
dZSTJ), (l^TUmi . and MUl to get 

ff 2 = Jj^i (4.47) 

and hence define ff Q = </ft-±ff as the invariant driving force for normal motion of the interface 
(superscript + is used for defining ff 7 and — for £ s ). 

To obtain a necessary condition for (14.43p . let A = (K^)" 1 (therefore A = (IfCj)" 1 ) locally 
at the point at which (|4Tl3]) is evaluated. With this choice for A, relations (f2^3|) . (f2l)5j) . (fOT) . 
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d%S!fl , and (f^TOTT) yield 

K 7 = I 7 , K 5 2 = Mi, (4.48) 

K 7 = (K 7 )- 1 ^ 7 , = (K 7 )" 1 ^ (4.49) 

ff 2 = J A .+ fi, f/ 2 = ^1^1, and (4.50) 

L 2 = (/i 7 )- 1 ^ 7 )- 1 ^^ 7 , (4.51) 

whereas HP, H a , and N a remain invariant. The above argument can be repeated with — superscript 
in place of + (and 5 in place of 7). We therefore obtain the necessary and sufficient condition for 
/ to be invariant under compatible changes in the reference configuration as 

/ = f(F*,M a ,M,W a ,(Ki)- 1 K a ,V*,N a ,U a ). (4.52) 

The proof of sufficiency is straightforward and therefore omitted. 

The list of arguments can be reduced on noting the definition (I3.3ip of M in addition to 

N 7 = m^MPN" 5 , f 5 = ^-P, f/ 7 = ^m~ l U 8 , (4.53) 

Jk+ Jk+ 

(K 7 )- 1 ^ 5 = 1 7 M + (K 7 ) _1 ]K 7 M, and (4.54) 

L «5 = ^Z m -iM" 1 L 7 M, (4.55) 
Jk+ 

where m = |M*| = j S (p)~^- The representation (I4.52h can then be equivalently written in the 
form 

/ = /(ff a , H 7 , M, H a , M, (K 7 ) _1 lk 7 , L 7 , N 7 , t/ 7 ). (4.56) 

This provides us with a complete and mutually independent set of variables that can be used as 
arguments for constitutive functions which are invariant under compatible transformation of the 
reference configuration. 

Under superimposed rigid body motions, all except H 7 and H 7 in the arguments of / remain 
invariant. They transform to QH 7 and QH 7 + QH 7 respectively, where Q G Orth + is arbitrary, 
cf. Subsection 12.51 For / to be invariant under superimposed rigid body motions we require 

/(H 7 ,H 7 , • • • ) = /(QH 7 , QH 7 + QH 7 , • • • ), (4.57) 

where dependence on other variables is suppressed. Choose Q = (M 7 ) T + N 7 (g>ii where R 7 is defined 
in K26h . With this choice QH 7 = U 7 and QH 7 + QH 7 = U 7 . The second of these follows on using 
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(M 7 ) T IR 7 = I 7 , (M 7 ) T ii = 0, U 7 N 7 = 0, and their normal time derivatives. Recalling C 7 = (IF) 2 
we thus have the necessary and sufficient condition for invariance under rigid body motions: 

/ = /(f ^ 7 ,!^,!: 7 ,!^,^ 7 )- 1 ^ 7 ,^^ 7 ,^ 7 ). (4.58) 

The proof for sufficiency is straightforward and therefore omitted. 

Assume for illustrative purposes that the conditions required by inverse function theorem are 
satisfied such that we can invert / = 0, where / is as given in f)4.58j) . to obtain 

W = ^(ff^C^M^M^K 7 )- 1 ^ 7 ^ 7 ^ 7 ), (4.59) 
M = ^(f a ,C 7 ,M,C 7 ,(K 7 )~ 1 ]K 7 ,L 7 ,N 7 ,C/ 7 ), and (4.60) 
(IK 7 )- 1 ^ 7 = S(ff Q ,C 7 ,M,C 7 ,M,L 7 ,N 7 ,£/ 7 ). (4.61) 

The following restrictions imposed by the material symmetry on f/ 7 , 1Z, and S ensue on recalling 
the pertinent discussion from Subsection 12.51 

jji = ^(ff^C^M^M, (IK 7 )" 1 ^ 7 ,^ 7 ,^ 7 ), (4.62) 
(G 7 ) T 71G 5 = TZ{f a , C 7 , M, C 7 , (IK 7 )- 1 ^ 7 , L 7 , N 7 , t/ 7 ), (4.63) 
and (G 7 ) T SG 7 = S{f a , C 7 , M, C 7 , M, L 7 , N 7 , C/ 7 ), (4.64) 

where (recall G a = t a G a , . and (l2~T03D ) 

C 7 = (G 7 ) T C 7 G 7 , M = (G 7 ) r MG 5 , (4.65) 
C 7 = (G 7 ) T C 7 G 7 , M = (G 7 ) T MG' 5 , (t 7 )"^ 7 = (G 7 ) T (IfC 7 )~ 1 lK 7 G 7 , (4.66) 
L 7 = (G 7 ) T L 7 G 7 , and N 7 = (G 7 ) T N 7 . (4.67) 

Here G 7 € £ 7 and G 5 G G s are the symmetry maps at the interface. The scalars f a and U a remain 
invariant since Jc a = 1 and |(G Q )*| = 1. 

Remark 4.5. If / is independent of N 7 and N* 5 then its dependence on C 7 , M, and K7 is only 
through 1 7 C 7 1 7 , 1 7 M1 <5 , and IK 7 1 7 , respectively. These follow immediately on taking the normal 
time derivative of C 7 N 7 = 0, M T N 7 = 0, MN 5 = 0, and K 7 N 7 = 0. 

Remark 4.6. A boundary-initial-value problem for %(X,t) and K(X, i) is specified through the 
coupled system of nonlinear partial differential equations given by (|2.53[) (with P substituted from 
(|4.5p ) and a flow rule of the type (cf. equation (102) of [To] ) 

K-^K = H(C H ,C H , a ), (4.68) 
47 



where Ch = H T H and a is the bulk dislocation density defined in (|3.22[) . The boundary data on 
S r are furnished by (|2.54p combined with the interfacial kinetic laws (|4.59p — ( |4.6ip . The problem 
is completed by supplying appropriate boundary conditions on dn r and initial conditions for x an d 
K at some fixed time. 

Anurag Gupta acknowledges the support of the initiation grant from Indian Institute of Technology, 
Kanpur, India. 



References 

[1] K. E. Aifantis and J. R. Willis. The role of interfaces in enhancing the yield strength of 
composites and polycrystals. Journal of Mechanics and Physics of Solids, 53:1047-1070, 2005. 

[2] B. A. Bilby. Types of dislocation source. In Report of Bristol Conference on Defects in 
Crystalline Solids, pages 124-133. London: Physical Society, 1955. 

[3] B. A. Bilby, R. M. Bullough, and D. K. de Grinberg. General theory of surface dislocations. 
Discussions of the Faraday Society, 38:61-68, 1964. 

[4] B. A. Bilby, R. M. Bullough, and E. Smith. Continuous distributions of dislocations: A new 
application of the methods of non-Riemannian geometry. Proceedings of the Royal Society at 
London A, 231:263-273, 1955. 

[5] R. M. Bullough and B. A. Bilby. Continuous distributions of dislocations: Surface dislocations 
and the crystallography of martensitic transformations. Proceedings of the Physical Society B, 
69:1276-1286, 1956. 

[6] P. Cermelli and M. E. Gurtin. The dynamics of solid-solid phase transitions. 2. Incoherent 
interfaces. Archive of Rational Mechanics and Analysis, 127:41-99, 1994. 

[7] P. Cermelli and M. E. Gurtin. On the kinematics of incoherent phase transitions. Acta 
Metallurgica et Materialia, 42:3349-3359, 1994. 

[8] P. Cermelli and M. E. Gurtin. On the characterization of geometrically necessary dislocations 
in finite plasticity. Journal of Mechanics and Physics of Solids, 49:1539-1568, 2001. 



48 



[9] P. Cermelli and S. Sellers. Multi-phase equilibrium of crystalline solids. Journal of Mechanics 
and Physics of Solids, 48:765-796, 2000. 

[10] J. W. Christian. The Theory of Transformations in Metals and Alloys. Pergamon (Elsevier), 
2002. 

[11] C. Davini. Elastic invariants in crystal theory. In J. M. Ball, editor, Material Instabilities in 
Continuum Mechanics and Related Mathematical Problems, pages 85-105. Oxford, 1988. 

[12] M. Epstein and M. Elzanowski. Material Inhomogeneties and their Evolution. Springer, 2007. 

[13] M. Epstein and G. A. Maugin. The energy-momentum tensor and material uniformity in finite 
elasticity. Acta Mechanica, 83:127-133, 1990. 

[14] N. A. Fleck and J. R. Willis. A mathematical basis for strain-gradient plasticity theory. Part 
II: Tensorial plastic multiplier. Journal of Mechanics and Physics of Solids, 57:1045-1057, 
2009. 

[15] P. Gudmundson. A unified treatment of strain gradient plasticity. Journal of Mechanics and 
Physics of Solids, 52:1379-1406, 2004. 

[16] A. Gupta, D. J. Steigmann, and J. Stolken. On the evolution of plasticity and incompatibility. 
Mathematics and Mechanics of Solids, 12:583-610, 2007. 

[17] A. Gupta, D. J. Steigmann, and J. Stolken. Aspects of the phenomenological theory of elastic- 
plastic deformation. Journal of Elasticity, 104:249-266, 2011. 

[18] M. E. Gurtin. The dynamics of solid-solid phase transitions. 1. Coherent interfaces. Archive 
of Rational Mechanics and Analysis, 123:305-335, 1993. 

[19] M. E. Gurtin. Thermomechanics of evolving phase boundaries in the plane. Oxford University 
Press, 1993. 

[20] M. E. Gurtin. Configurational forces as basic concepts of continuum physics. Springer, 2000. 

[21] M. E. Gurtin and A. I. Murdoch. A continuum theory of elastic material surfaces. Archive of 
Rational Mechanics and Analysis, 57:291-323, 1975. 



49 



[22] M. E. Gurtin and A. Struthers. Multiphase thermomechanics with interfacial structure. 3. 
Evolving phase boundaries in the presence of bulk deformation. Archive of Rational Mechanics 
and Analysis, 112:97-160, 1990. 

[23] M. E. Gurtin, A. Struthers, and W. O. Williams. A transport theorem for moving interfaces. 
Quarterly of Applied Mathematics, 47:773-777, 1989. 

[24] C. Herring. Some theorems on the free energies of crystal surfaces. Physical Review, 82:87-93, 
1951. 

[25] F. J. Humphreys and M. Hatherly. Recrystallization and related annealing phenomena, Second 
edition. Pergamon (Elsevier), 2004. 

[26] W. C. Johnson and J. I. D. Alexander. Interfacial conditions for thermomechanical equilibrium 
in two-phase crystals. Journal of Applied Physics, 59:2735-2746, 1986. 

[27] C. C. Koch, I. A. Ovidko, S. Seal, and S. Veprek. Structural Nanocrystalline Materials: 
Fundamentals and Applications. Cambridge University Press, 2007. 

[28] K. Kondo. Non-Riemannian geometry of imperfect crystals from a macroscopic viewpoint. 
In K. Kondo, editor, Memoirs of the Unifying Study of the Basic Problems in Engineering 
and Physical Science by Means of Geometry, volume I, pages 458-469. Gakujutsu Bunken 
Fukyu-Kai, Tokyo, 1955. 

[29] K. Kondo. On the analytical and physical foundations of the theory of dislocations and yielding 
by the differential geometry of continua. International Journal of Engineering Science, 2:219- 
251, 1964. 

[30] E. Kroner. Continuum theory of defects. In R. Balian, M. Kleman, and J. Poirier, editors, Les 
Houches, Session XXXV, 1980 - Physique des defauts, pages 215-315. North-Holland, New 
York, 1981. 

[31] F. C. Larche and J. W. Cahn. Thermo chemical equilibrium of multiphase solids under stress. 
Acta Metallurgica, 26:1579-1589, 1978. 

[32] F. C. Larche and J. W. Cahn. The interaction of composition and stress in crystalline solids. 
Acta Metallurgica, 33:331-357, 1985. 



[33] P. H. Leo and R. F. Sekerka. The effect of surface stress on crystal-melt and crystal-crystal 
equilibrium. Acta Metallurgies 37:3119-3138, 1989. 

[34] M. A. Meyers. Dynamic Behavior of Materials. John Wiley & Sons Inc., 1994. 

[35] F. R. N. Nabarro and M. S. Duesbery, editors. Dislocations in Solids, volume 11. Elsevier, 
2002. 

[36] W. Noll. Materially uniform simple bodies with inhomogeneities. Archive of Rational Mechan- 
ics and Analysis, 27:1-32, 1967. 

[37] R. Penrose. A generalized inverse for matrices. Mathematical Proceedings of the Cambridge 
Philosophical Society, 51:406-413, 1954. 

[38] M. Silhavy. The Mechanics and Thermodynamics of Continuous Media. Springer, 1997. 

[39] M. Silhavy. Phase transitions with interfacial energy: Convexity conditions and the existence 
of minimizers. In J. Schroder and P. Neff, editors, CISM Course on Applications of Poly-, 
Quasi-, and Rank-One Convexity in Applied Mechanics, volume 516, pages 177-240. Springer, 
Wien and New York, 2010. 

[40] M. Silhavy. Phase transitions with interfacial energy: Interface null Lagrangians, polycon- 
vexity, and existence. In K. Hackl, editor, IUTAM Symposium on Variational Concepts with 
Applications to the Mechanics of Materials, pages 233-244. Springer, Dordrecht, 2010. 

[41] M. Silhavy. Equilibrium of phases with interfacial energy: A variational approach. Journal of 
Elasticity, 105:271-303, 2011. 

[42] D. J. Steigmann and R. W. Ogden. Elastic surface-substrate interactions. Proceedings of the 
Royal Society at London A, 455:437-474, 1999. 

[43] A. P. Sutton and R. W. Balluffi. Interfaces in Crystalline Materials. Oxford University Press, 
2009. 

[44] T. Y. Thomas. Plastic Flow and Fracture in Solids. Academic Press, 1961. 

[45] C. Truesdell and R. A. Toupin. The classical field theories. In S. Fliigge, editor, Handbuch der 
Physik, volume III/l, pages 226-793. Springer- Verlag, 1960. 



51 



